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High-ion absorption in seven GRB host galaxies at z=2-4* 
Evidence for both circumburst plasma and outflowing interstellar gas 
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ABSTRACT 

Aims. We use VLT/UVES high-resolution optical spectroscopy of seven GRB afterglows (at zgrb between 2.20 and 
3.97) to investigate circumburst and interstellar plasma in the host galaxies of the bursts. 

Methods. GRBs 050730, 050820, 050922C, 060607, 071031, and 080310 were each detected by the Swift satellite. With 
the minimum possible time delay (as short as eight minutes), follow-up optical spectroscopy at 6.0 kms -1 resolution 
began with UVES. We present Voigt-profile component fits and analysis of the high-ion absorption detected in these 
spectra and in the pre-Swtft UVES spectrum of GRB 021004, at velocities between zgrb-5 000 kms -1 and zgrb- 
Results. We identify several distinct categories of high- ion absorption at velocities close to zgrb: (i) Strong high-ion 
absorption at zgrb itself is always seen in Ovi, Civ, and Siiv, usually (in six of seven cases) in Nv, and occasionally 
in Siv and S VI. Three of the cases show logiV(Nv)>14 in a single-component, suggesting a circumburst origin, but we 
cannot rule out an interstellar origin. Indeed, using the non-detection of S IV* at zgrb toward GRB 050730 together 
with a UV photo-excitation model, we place a lower limit of 400 pc on the distance of the S iv-bearing gas from the 
GRB. (ii) Complex, multi-component Civ and Siiv profiles extending over 100-400 kms -1 around zgrb are observed 
in each spectrum; these velocity fields are similar to those measured in C IV in damped Lyman-a systems at similar 
redshifts, suggesting a galactic origin, (iii) Asymmetric, blueshifted, absorption-line wings covering 65-140 kms" 1 are 
seen in the Civ, Siiv, and Ovi profiles in four of the seven GRB afterglow spectra. The wing kinematics together 
with the "Galactic" Civ/Siiv ratios measured in two cases suggest that the wings trace outflowing interstellar gas 
in the GRB host galaxies, (iv) High-velocity (HV; 500-5 000 kms -1 relative to zgrb) components are detected in six 
of the seven spectra; these components are not necessarily a single homogeneous population, and many of them may 
arise in unrelated foreground galaxies. However, in the cases of GRBs 071031 and 080310, the ionization properties of 
the HV components (very high Civ/Siiv ratios and absence of neutral-phase absorption in Sill or Cn) are suggestive 
of a circumburst origin; models of Wolf-Rayet winds from the GRB progenitors can explain both the kinematics and 
ionization level of these HV components. 

Key words, gamma rays: bursts - galaxies: halos - galaxies: high-redshift - galaxies: ISM - stars: Wolf-Rayet 



1. Introduction 

For many years the only backlights available for study- 
ing the distant Universe through absorption-line spec- 
troscopy were AGN. Since the d iscovery that gamma-ray 
bursts (GRBs) are e xtragalactic (|van Paradiis et aLlll997t 
iMetzger et al.1 [l997h . GRBs have also been employed as 
backlights, with the advantages of extremely high lumi- 
nosity and power-law continua, but the disadvantage of 
rapid temporal fading, implying that rapid-response obser- 
vations are necessary for their study. Not only do anal- 
yses of GRB afterglow spectra allow intervening systems 
to be studied, they also enable the properties of the in- 
terstellar gas in the host galaxy of the GRB to be inves- 
tigated, particularly its kinematics, chemical abundances, 
dust content, molecular gas content, and ionized gas con- 



tent (ICastro et alj2003tlKlose et al.ll2004tlFiore et alj2005t 



i Savagfio et all 120031 . ISavaglio fc Fall 120041 
2006: iFvnbo et all 1200a " 

Many GRB spectra show 
sorption at zgrb> i-e 



Pcnp rase et ahl 
Prochaska et all " 120061 l2007aD . 
damped Lyman-a (DLA) ab- 
t hey show a neutra l gas column 



densi t y log iV(Hi)>20. 3 (iJensen et al.l l200lt Irliorth et al.l 



2003; Jakobsson et al 



2006; IVreeswiik et all 2004 



Starling et al.ll2005HBerger et alj|2006t I Watson et al.l l2006 



iProchaska et al. l l2007bD . typically with higher Hi column 
densities than ar e seen in DLAs toward QSOs (QSO- 
DLAs, reviewed bv lWolfe et aHl2005l ). The higher JV(H i) in 
GRB-DL As may be re l ated to their intense star-formation 
activity (jBloom et alj I2002D . or to the GRB sight-lines 
passing prefer entially through the inner regions of th e 
host galaxies (jFvnbo et al.l 120081 ; IProchaska et aTll2008aD . 
Photometric studies of GRB host g alaxies have shown them 



* Based on observations taken under Programme IDs 
070.A-0599, 070.D-0523, 075.A-0385, 075.A-0603, 077.D-0661, 
and 080.D-0526 with the Ultraviolet and Visual Echelle 
Spectrograph (UVES) on the Very Large Telescope (VLT) Unit 
2 (Kueyen) at Paranal, Chile, operated by ESO. 



IWiersema et~a!]|2007t IThone et al l 120071 l2008al). with 

(|Sa\ 



no 



to be star-forming, dwarf galaxies (IChristensen et al] [2004; 

evidence that the hosts are peculiar (|Savaglio et al.l l2008) 

In this paper we use rapid-response mode spectra of 
seven high-redshift long-duration GRB afterglows to inves- 
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tigate high-ion absorption (in Ovi, Nv, Civ, Siiv, Siv, 
and Svi) in the ISM of the host galaxy, as well as in the 
circumburst medium immediately surrounding the GRB. 
Observations of ionized interstellar gas in galactic envi- 
ronments provide important constraints on several phys- 
ical processes, including star formation and subsequent 
feedback, galactic winds, interactions with satellite galax- 
ies, and accretion. In the pve-Swift era, studying high 
ions in galaxy environments at high redshifts was compli- 
cated by the difficulty of associating individual absorbers 
with individual galaxies. There are several indirect ways 
to explore high-ion absorption in high-rcdshift galaxies. 
Firstly, outflowing Civ and Siiv is seen in the compos- 
ite spectrum of high -redshift Lyman-break galaxies (LBGs; 
IShaplev et al.|[2003h . Secondly, there is a statistical corre- 
lation between the redshifts of C iv absorbers seen in QSO 
spectra and the redshifts of ga laxies near the QSO lines- 
of-sight (jAdelberger et al.ll2005l ). And thirdly, high- ion ab- 
sorption is detec te d in all DLAs (TWolfe fc Procha ska 2000; 
iFox et al.|[2007aUb1 : iLehner et alj |2008f). which trace galac- 
tic structures. However, because of the faintness of high- 
redshift gal axies, only in a sing l e pub lished case (the lensed 
LBG cB58; lPettini et al.|[2000l . I2002D are interstellar high- 
ion absorption lines seen at moderate-to-high resolution in 
the spectra of individual galaxies at z>2. 

Now, with the ability to be on-target within a few min- 
utes of a Swift GRB trigger, a new means to directly probe 
the interstellar (as well as circumburst) medium in individ- 
ual high-redshift g alaxies is available (see lFiore et al.l[2005t 
iD'Elia et al1l2007t IChen et alj|2007h . High-resolution opti- 
cal spectroscopy of GRB afterglows now routinely produces 
spectra with high-enough signal-to-noise to uncover the ve- 
locity sub-structure in the absorbing gas. Evidence has even 
been found for time- variation in the strength of absorption 
in fine-structure lines at zg rb (|Dessauges-Zavadskv et al] 
2006; IVreeswiik et al.|[2007l ). In this paper, we present a 
survey of high-ion absorption in the highest-quality GRB 
optical afterglow spectra taken to date. We pay particular 
attention to the high-ion kinematics, since velocity mea- 
surements can constrain the presence of gaseous outflows, 
either from the GRB progenitor or from the host galaxy. 
We structure this paper as follows. In §2 we discuss the 
GRB observations, data processing, and our line profile fit- 
ting procedures. In §3 we describe the observational prop- 
erties of the high-ion absorption at velocities near zgrb 
in each afterglow spectrum. We discuss four distinct ob- 
servational categories of high-ion absorption in the GRB 
afterglow spectra in §4, and we present a summary in §5. 

2. GRB observations 

We formed our sample of seven GRB afterglow spectra by 
selecting all GRBs at z > 2 with follow-up high-resolution 
UVES3 optical spectroscopy available as of March 2008. 
Lower redshift GRBs were not considered since no infor- 
mation on the high ions (particularly O vi) is available in 
these cases. With the exception of GRB 021004, which was 
detected by the High-Energy Transient Explorer (HETE- 
2) satellite (jRicker et al]|2003l ). each of the GRBs in our 
sample was detected by the Burst Alert Teles cope ( BAT) 
on-board NASA's Swift satellite (jGehrels et al.ll2004l ). The 
seven GRBs are 021004, 050730, 050820, 050922C, 060607, 

1 UVES is described in lDekker et~aTl (|2000Q . 



071031, and 080310, named according to the UT date of 
their detection (yymmdd). Following each trigger, Swift be- 
gan follow-up observations with the XRT and UVOT in- 
struments, to determine an error circle of w2.5" around 
the source. 

After the Swift localization was distributed through the 
GCN (GRB Coordinates Network), rapid-response mode 
(RRM) UVES observations began on each GRB afterglow 
with the minimum possible time delay, provided the target 
was observable from Paranal. Three of the seven GRB af- 
terglows in our sample were observed in this mode. RRM 
observations involve the immediate interruption of any cur- 
rent UVES exposure in order to slew to and acquire the 
new target. No human intervention is required, except for 
alignment of the target on the spectrograph slit. On two 
occasions, a sub-ten minute response was achieved. These 
observations demonstrate the power of the RRM concept, 
and the success of its implementation on ESO telescopes. 
The four afterglow spectra that were not taken in the RRM 
mode are GRBs 021004, 050820, 050730 and 050922C, 
which either occurred before the RRM was implemented 
or were unobservable from Paranal at the time of the trig- 
ger. In each of these cases, the UVES observations began 
at the earliest possible time in the form of a traditional 
Target-of-Opportunity (ToO) observation. A log of the ob- 
servations, together with a summary of the basic properties 
of the GRBs, is given in Table 1. 

The spectra were taken with 2x2 binning and a 1.0" 
slit, providing R w 50 000 (FWHM=6.0 krns" 1 ) spectra 
in «2 kms" 1 pixels0 The UVES data were reduced with 
a customized version of the MIDAS reduction pipeline 
()Ballester et al] [2000) , with the wavelength scale subse- 
quently corrected into the vacuum heliocentric frame. For 
the latest GRB in our sample (GRB 080310), we also re- 
duced the data independently with version 3.4.5 of the CPL 
(Common Pipeline Library) pipelinfl A detailed compar- 
ison of line profiles found the results from the two pipeline 
versions to be essentially identical. Individual exposures 
were combined with a least-squares weighting to produce 
a combined spectrum. In the case of GRB 071031, we also 
maintain and analyze the individual exposures to look for 
time-series variation (see §3.6.1). Continua were fit locally 
to each absorption line of interest. Before displaying and 
fitting the Ovi, Cm, Nv, and Svi profiles, which fall in 
low S/N regions of the spectrum, we perform a further re- 
binning, generally by five pixels, but occasionally by other 
factors, as described in the captions to Figures 1 to 4. These 
differing rcbinning factors are driven by the desire to max- 
imize the signal-to-noise ratio and to investigate the signif- 
icance of the high-ion detections, even if some resolution is 
compromised. 

We searched in each spectrum for high-ion ab- 
sorption over the velocity range zgrb~5 000 kms -1 to 
zgrb+2000 kms -1 , where zgrb is defined by the posi- 
tion of strongest absorption in the low- ionization lines. This 
wide range was chosen because of the possibility of detect- 
ing high-velocity outflow features driven by the burst pro- 
genitor. We make a distinction between low- velocity (LV) 
absorption components seen within several hundred kms -1 



2 Though the nominal spectral resolution of UVES in this 
mode is 43 000, we find that a higher resolution of f»50 000 is 
achieved in practice, due to variations in the seeing conditions. 

3 See manual at [ftp://ftp.eso. org/ pub/dfs/pipeli nes/uves/| 
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UT(UVES)^ 
(hh:mm) 



GCN' 
Refs" 



GRB 
(yymmdd) 



UT(Trig^ 



8t° 
(hh:mm) 



ESQ ID 



00 



2GRB 10gJV H ^ [Z/H] e 



Science^ 
Refs. ~ 



021004 


12:06:13 


01:37 


13:31 


7200 


070.A-0599.2 


2.3290 


19.00 


0.0 


1,2 


3,4,5,6,7,8 


050730 


19:58:23 


00:07 


04:09 


6000 


075.A-0603" 


3.9686 


22.15 


-2.3 


9,10 


8,11,12,13,14,15,16 


050820 


06:34:53 


07:08 


00:34 


6043 


075.A-0385 


2.6147 


21.00 


-0.6 


17,18 


8,13,15,16 


050922C 


19:55:50 


23:42 


03:47 


6000 


075.A-0603 


2.1990 


21.55 


-2.0 


19,20 


8,16,21 


060607 


05:12:13 


05:20 


00:08 


11980 


077.D-0661 


3.0749 


16.80 


0.0 


22,23 


8 


071031 


01:06:36 


01:15^ 


00:09 


9480 


080.D-0526 


2.6922 


22.15 


-1.7 


24,25 




080310 


08:37:58 


08:5~ 


00:13 


4680 


080.D-0526 


2.4274 


18.80 


-1.4 


26,27 





a UT of trigger by the BAT instrument on-board Swift. Exception: GRB 020104, detected by WXM on-board HETE-2. 

b UT of start of first UVES exposure (after target acquisition). 

c Time delay between satellite trigger and start of UVES exposure. 

d Total UVES exposure time over all setups. 

E Hi column density and metallicity of GRB-DLA taken from reference 8, except 071031 and 080310 values, from this paper. 
' We list here the GCN first reporting the burst, and the GCN first reporting the UVES spectrum. 
9 We list references relevant to the study of the host galaxy ISM. 
h Also 070.D-0523. 

1 Due to a difficult target acquisition, the first fully aligned exposure of GRB 071031 began at 01:25. 
REFE R ENCES: llShirasaki et all d2002ft; 2lSavaglio et all J2002D: 3lMe>ller et all (l2002f>: 4lMirabal et all (120031): 5lSchaefer et all 
d2003fl: 6lFiore et all (120051 1: 7 lLazzati et all (120061 1: 8 IProchaska et all d2008bf>: 9lHolland et all J2005fi: lOlD'Elia et all (l2005al T: 
11 IStarling et alj d2005lj~ 12 IChen et all (1200511 : 13 IChen et all (1200711 :! Tl D'Elia et all 11200711: 1 5 IProchaska et al l il2007all~16 



Prochaska ct al 



.2007bl1 : 1 7 iPaee et all (120051): 1 8 ILedoux et all (1200511: 19 iNorris et all (1200511: 20 iD'Elia et"aiT (120051: 
200611 : 22 IZiaeepour et all ll2006Tl : 23 ILedoux et all (|2006t 1: 24 IStroh et"all (|2007l 1: 25 | Ledoux et all (12007 



Jakobsson ct al. , _ 

Cummings et all (|2008l ); 27 lVreeswiik et all (|2008Q 
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of zgrb, which are seen in all cases, and high- velocity (HV) 
components at 500-5000 kms" 1 from zqrb, which are seen 
in 6/7 sight lines, and which are clearly separated from the 
LV absorption. When absorption was found, we used the 
VPFIT software package^ to fit the high-ion components 
in each GRB spectrum with a series of Voigt profiles. Each 
high ion was fit independently, even species with similar 
ionization potentials. The VPFIT software accounts for in- 
strumental resolution, and returns the redshift, line width, 
and column density (z, b, and logiV) for each component 
together with the associated errors. We convert the compo- 
nent redshift to velocity, using zqrb as the zero point. By 
summing over the component column densities, and adding 
the component errors in quadrature, we form the total col- 
umn density and its error in each absorber. In cases where 
the high-ion absorption is saturated, the errors on the total 
column density are large. The fit solutions are non-unique, 
since the number of components to be fit has to be speci- 
fied manually. This number is often larger than ten, due to 
the complexity of the absorption-line profiles. Nonetheless, 
our fitting technique has been tested by comparing the to- 
tal model column densities with th e values from dire ct ap- 
parent optical depth integrations (|Fox et all l2007al ). The 
successful outcome of these tests supports the reliability of 
the fit results. The details of the fits are given in Tables 2 
to 8, for GRBs 021004, 050730, 050820, 050922C, 060607, 
071031, and 080310, respectively. Ato mic data (rest wave- 
lengths and /-values) were taken from iMortonl ((2003). 



Table 2. GRB 021004 z=2.3290: Voigt profile fits 



3. Description of the high-ion absorption 

High-ion absorption within several hundred kms -1 of zqrb 
is found in all seven GRB spectra in our sample. The de- 
tected lines are Civ and Siiv (seen in 7/7 cases), Nv 



Ion 


«o 

(kms -1 ) 


b 

(kms -1 ) 


log TV 
(N in cm" 2 ) 


log Af (total) 
(TV in cm -2 ) 


Ovi 


-25±8 


18±10 


14.46±0.22 


14.95±0.20 




26±7 


22±8 


14.78±0.22 




Nv 


-2903±2 


19±2 


14.03±0.05 


14.03±0.05 


Nv 


-10±2 


28±3 


14.53±0.06 


14.53±0.06 


Civ 


-2966±4 


12±6 


12.97±0.19 


14.73±0.08 




-2903±2 


30±2 


14.59±0.04 






-2720±2 


8±2 


13.36±0.09 






-2670±2 


31±3 


14.03±0.04 






-2599±8 


24±11 


13.10±0.18 




Si iv 


-2903±3 






ft;14.82* 




-2882±45 


58±40 


13.53±0.42 






-2659±4 


20±5 


13.20±0.10 




Siv 


0±4 


45±5 


15.25±0.20 


15.25±0.20 



Available at http://www.ast.cam.ac.uk/~rfc/vpfit.html 



NOTE: the symbol * denotes a saturated component; no re- 
liable column density or line width can be extracted for the 
component, and the total column density in that ion is un- 
certain. NOTE: Total column densities (summed over compo- 
nents) are reported separately for the HV and LV absorption. 
NOTE: No Civ or Siiv fit was attempted to the absorption 
within 1000 kms -1 of zgrb, due to complete saturation. 



(6/7), Ovi (6/7), Siv (5/7), and Svi (2/7). The high- 
ion absorption-line profiles in the velocity range within 
200 kms -1 of zqrb are shown for each of our seven GRB 
afterglow spectra in Figures 1 to 4. In addition to the high 
ions, we include on each figure absorption lines that trace 
the neutral gas, typically Sin A1260.422, Cn A1334.532, or 
Fen A1608.451. 

In cases where we do not report a detection for Ovi, 
Siv, and Svi (e.g. Ovi toward GRB 060607), it is because 
the profiles are blended, not because there are actual non- 
detections in these ions. The only clear non-detection at 
•^GRB in any of the high ions mentioned above is in Nv 
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Table 3. GRB 050730 z=3.9686: Voigt profile fits 



GRB 021004 z=2.3290 log N(H I)=19.00 



Ion 


Vo 

(kms" 1 ) 


b 

(kms -1 ) 


log TV 
(TV in cm -2 ) 


log JV (total) 
(TV in cm -2 ) 


Ovi 


-88±5 


34±5 


14.52±0.08 


«15.94* 




-54±2 


* 








1±8 


36±6 


14.97±0.15 






31±6 


* 


* 




Nv 


- 1±2 


7±3 


13.54±0.15 


13.80±0.14 




26±5 


14±8 


13.46±0.18 




Civ 


-1577±2 


4±2 


13.38±0.22 


14.52±0.18 




-1567±2 


5±4 


13.43±0.23 






-1553±2 


5±2 


13.66±0.17 






-1536±2 


12±2 


13.97±0.10 






-1516±2 


34±34 


13.81±0.37 






-1480±2 


11±3 


13.86±0.17 




Civ 


-168±2 


14±3 


13.12±0.06 


«15.24* 




-119±18 


* 


* 






-98±7 


* 


* 






-58±2 


17±3 


14.81±0.10 






4±2 


21±2 


15.01±0.06 




Si iv 


-1575±2 


4±2 


12.57±0.06 


13.94±0.10 




-1553±2 


4±2 


13.03±0.05 






-1542±3 


26±3 


13.38±0.08 






-1529±2 


7±4 


12.45±0.38 






-1499±2 


9±4 


12.61±0.21 






-1472±2 


12±2 


13.50±0.03 






-1442±2 


13±2 


13.02±0.04 




Si iv 


-104±6 


* 


* 


«14.35* 




-63±50 


* 


* 






-35±50 




* 






1±2 


10±2 


13.94±0.12 




Siv 


-53±2 


28±4 


14.51±0.06 


«15.31* 




-55±2 


* 








1±2 


8±2 


14.65±0.06 






21±2 


4±3 


13.57±0.15 




Svi 


- 3±2 


17±2 


14.16±0.06 


14.16±0.06 



o 



NOTE: the symbol * denotes a saturated component; no reli- 
able column density or line width can be extracted for the com- 
ponent, and the total column density in that ion is uncertain. 
NOTE: Total column densities (summed over components) are 
reported separately for the HV and LV absorption. 



toward GRB 060607. In two cases our data cover the wave- 
length of Nevm AA770.409, 780.324 at z G rb, but the dou- 
blet lies below the Lyman limit where no flux is transmit- 
ted. Pv AA1117.977, 1128.008 is not detected in any of 
the seven spectra. Cm A977.020 is detected in one case 
(GRB 071031), but in the other sight lines the S/N is too 
low at the observed wavelength (far into the blue) for the 
C in data to be useful. Absorption in S IV* A1072.973 is not 
observed in any of the seven spectra. Finally, with the ex- 
ception of GRB 060607 with its unusually low TV(H i), Si in 
A1206.500 lies in the strong damping wings of Lyman-a, so 
no useful information on this ion at zqrb is available. 

The detection rate of Nv in our samp le is identical to 
that reported bv iProchaska et all (|2008bl ). who find Nv 
at £grb in 6/7 GRB afterglow spectra, but since five sys- 
tems are common to both samples, this is not surprising. 
N v has also been det ected at zgrb in th e afterglow spec- 
trum of GRB 060206 (jThone et al.ll2008bf h Other than the 
O vi absorption in the GRB 050730 s pectrum, which has 
been published bv lD'Elia et all (|2007l ). no reports of Ovi 
at zqrb in GRB afterglow spectra exist, and no previous 
reports of S iv or S vi exist either. Indeed, very few re- 
ports of interstellar S vi absorption exist anywhere in the 




-200 
v (km s~') 



200 



Fig. 1. Normalized high- ion and low-ion absorption- line 
profiles on a velocity scale relative to zqrb for GRB 021004. 
Voigt profile fits are shown in blue (solid line), with the cen- 
ter of each component in the fit marked with a red dash. 
Dark shading denotes the strong absorption components at 
zgrb- All data have been rebinned by four pixels (eight in 
the case of Ovi). We do not fit the Civ or Si IV doublets 
because of saturation. Absorption components at higher ve- 
locity relative to zqrb are shown in Figure 5. 
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GRB 050730 z=3.9686 log N(H I)=22.15 GRB 050820 z=2.6147 log N(H I)=21.00 




Fig. 2. Normalized high- ion and low- ion absorption-line profiles on a velocity scale relative to zgrb for GRBs 050730 and 
050820. Voigt profile fits are shown in blue (solid line), with the center of each component in the fit marked with a red 
dash. Dark shading denotes the regions of strong absorption at zgrb itself, and light shading denotes the negative- velocity 
wings (GRB 050730 only). The Ovi, Nv, and Svi data have been rebinned by five pixels. Absorption components at 
higher velocity relative to zgrb are shown in Figure 5. 
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Fig. 3. Normalized high-ion and low-ion absorption-line profiles on a velocity scale relative to zqrb for GRBs 050922C 
and 060607. Voigt profile fits are shown in blue (solid line), with the center of each component in the fit marked with a 
red dash. Dark shading denotes the regions of strong absorption, and light shading denotes the negative-velocity wing 
(GRB 060607 only). The Siv and Svi profiles are blended in each of these sight lines. The O vi and Nv data have been 
rebinned by eight pixels for GRB 050922C, and the Nv data by five pixels for GRB 060607. Absorption components at 
higher velocity relative to zgrb are shown in Figure 5. 
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Fig. 4. Normalized high-ion and low-ion absorption-line profiles on a velocity scale relative to zqrb for GRBs 071031 
and 080310. Voigt profile fits are shown in blue (solid line), with the center of each component in the fit marked with a 
red dash. Dark shading denotes the regions of strong absorption, and light shading denotes the negative-velocity wings. 
The Ovi and Nv data have been rebinned by four pixels for GRB 071031 and by five pixels for GRB 080310. For Svi, 
the rebin factors are five and ten, respectively. Absorption components at higher velocity relative to zqrb are shown in 
Figure 5. 
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Fig. 5. Seven high- velocity (500-5 000 kms -1 from 2grb) high-ion absorbers are seen in the seven GRB afterglow spectra 
in our sample; these HV absorbers are shown here (GRB 021004 shows two HV absorbers, GRB 050922C shows none, 
and the other five GRB spectra show one each). Voigt profile fits are only shown (with blue solid lines) for cases where we 
judge that reliable column densities can be extracted. The properties of these HV absorbers are diverse, and it is unlikely 
that they represent a single population. Several of these HV absorbers are seen in both neutral-phase and high-ion-phasc 
absorption; others (e.g. the HV components toward GRB 071031 and GRB 080310) are seen only in high-ion transitions. 
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Table 4. GRB 050820 2=2.6147: Voigt profile fits Table 5. GRB 050922C z=2.1990: Voigt profile fits 



Ion 




b 


log TV 


log JV (total) 


Ion 


«o 


b 


log TV 


logTV(total) 
(TV in cm" 2 ) 




(kms" 1 ) 


(kms" 1 ) 


(TV in cm -2 ) 


(TV in cm -2 ) 




(kms" 1 ) 


(kms" 1 ) 


(TV in cm" 2 ) 


Nv 


-109±2 


9±3 


13.17±0.08 


13.17±0.08 


Ovi 


65±2 


24±2 


14.44±0.05 


14.44±0.05 


Civ 


-3655±2 


10±2 


13.02±0.04 


13.13±0.06 


Nv 


63±2 


16±3 


13.76±0.07 


13.76±0.07 




-3631±2 


10±3 


12.50±0.11 




Civ 


-20±2 


16±2 


13.39±0.03 


15.69±0.35 


Civ 


-268±6 


65±8 


13.42±0.05 


14.91±0.11 




25±2 


10±2 


13.21±0.08 






-247±2 


9±2 


13.87±0.03 






63±2 


11±2 


15.67±0.30 






-223±2 


9±2 


13.53±0.03 






96±2 


8±2 


13.63±0.07 






-167±2 


9±2 


13.13±0.07 






140±2 


18±2 


13.73±0.03 






-146±2 


9±2 


13.58±0.04 






201±2 


10±2 


13.18±0.04 






-109±2 


18±2 


14.02±0.02 




Si iv 


27±2 


7±2 


13.22±0.04 


14.49±0.16 




-62±2 


15±2 


13.74±0.06 






62±2 


12±2 


14.40±0.11 






-15±2 


25±2 


14.55±0.04 






85±2 


5±3 


13.13±0.24 






84±2 


35±2 


14.06±0.02 






98±2 


6±2 


13.28±0.05 




Si iv 


-249±2 


5±2 


12.23±0.05 


14.35±0.18 




127±2 


4±2 


12.42±0.07 






-222±2 


8±3 


11.98±0.10 






141±2 


4±2 


12.70±0.05 






-170±2 


8±2 


12.72±0.04 






198±2 


7±2 


12.53±0.07 






-145±2 


10±2 


13.05±0.06 
















-122±3 
-57±2 


13±4 
7±2 


12.67±0.14 
12.74±0.14 




Table 6. GRB 060607 z=3.0749: Voigt profile fits 




-37±3 


46±3 


13.72±0.08 
















-12±2 


20±2 


14.00±0.04 




Ion 


«o 


b 


log TV 


logTV(total) 
(TV in cm" 2 ) 




64±2 


16±2 


12.96±0.06 






(kms" 1 ) 


(kms" 1 ) 


(TV in cm" 2 ) 




91±2 


11±3 


12.82±0.09 




Civ 


-1849±2 


9±2 


13.70±0.03 


13.83±0.05 




111±2 


4±3 


12.14±0.20 






-1824±2 


18±2 


13.23±0.05 






147±2 


26±3 


13.23±0.05 




Civ 


-88±3 


31±4 


13.00±0.05 


15.94±0.32 




168±2 


4±3 


12.05±0.19 






-45±2 


10±2 


13.03±0.06 






220±2 


4±2 


12.39±0.13 






- 5±2 


19±2 


14.26±0.02 






202±6 


17±12 


12.32±0.29 






19±2 


4±2 


15.93±0.28 




Siv 


-145±2 


25±3 


14.04±0.04 


14.91±0.07 


Si iv 


-1848±2 


11±2 


13.12±0.04 


13.28±0.06 




-67±2 


4±3 


13.30±0.10 






-1817±2 


12±3 


12.76±0.07 






- 5±2 


4±2 


14.04±0.12 




Si iv 


-47±2 


4±2 


12.31±0.12 


13.54±0.09 




- 8±2 


32±2 


14.76±0.03 






-30±3 


8±5 


12.48±0.23 




S vi 


-109±3 


20±4 


13.45±0.07 


14.09±0.08 




- 6±2 


12±2 


13.30±0.06 






-49±3 


14±3 


13.54±0.11 






16±2 


10±2 


12.98±0.07 






-11±3 


17±3 


13.77±0.07 




NOTE: Total column densities (summed over components) are 



NOTE: Total column densities (summed over components) are 
reported separately for the HV and LV absorption. NOTE: S VI 
is only detected in A933 (A944 is blended), so the detection is 
uncertain. Nonetheless the velocity structure corresponds closely 
to C iv. 



high-z (IChurchill fc Charlton! 1 1 999t iLevshakov et al.l l2003) 
or low-z ( Savage et al. 2005b) Universe. Interstellar Siv 



absorption-line detections have (to our knowledge) only 
been reported in qua sar-intrinsic systems (|Arav et al.lll999t 
iGangulv et all I2006D . Finally, our non-detections of ab- 
sorption in Siv* A1072.973 contrast with detections of 
this line in the win ds of OB-type stars (|Massa et aLl feOOS: 
iLehner et ai1l2003f) . 

The absorption features observed at and near zgrb can 
be classified into various types. First, strong absorption is 
seen exactly at zgrb in the form of saturated Ovi, Civ, 
and Si iv lines, usually with accompanying N v, and (where 
the data are unblended) also in S iv and S VI. Second, com- 
plex multi-component C IV and Si IV profiles with total ve- 
locity widths of 100-400 kms -1 are seen in all seven GRB 
afterglow spectra. Third, 4/7 GRB afterglow spectra show 
asymmetric negative- velocity absorption-line wings extend- 
ing over velocity regions of 65-140 kms -1 . These wings are 
seen clearest in C iv and Si iv, but are also present in O vi. 
Fourth, 6/7 cases show HV components in Civ, and oc- 
casionally also in neutral-phase lines such as Si n and C II. 
Whereas the first three of these categories are visible in 
Figures 1 to 4, we present a separate plot (Figure 5) show- 



reported separately for the HV and LV absorption. 



ing the HV high-ion absorbers detected in our sample. We 
return to a discussion of these categories of absorption in 
§4, after describing the details of the high-ion absorption in 
each spectrum. Unless otherwis e stated, solar abundances 
were taken from|LodderJ (|2003l ). 

3.1. GRB 021004 z GRB =2. 3290 

Extremely strong C iv absorption is seen at zgrb , with a 
fully saturated trough spreading over 300 kms -1 . Siiv is 
also present and less saturated than C IV, but we cannot dis- 
cern the true component structure in Si IV because of low 
S/N. A component centered at kms" 1 is seen in Ovi, 
Nv, and Siv, with absorption also detected in Aim, Sin, 
Cn, and Fen. Siv* A1072.973 is not detected at z G rb; we 
measure log TV(S iv*)<14.19 (3<r) over the velocity range 
-60 to 60 kms" 1 , and log TV(S iv)=15.30±0.29, so that 
TV(SlV*)/TV(Siv)<0.08. Two HV absorbers are observed: 
one centered at —2 900 kms" 1 relative to zqrb, the other 
centered at «— 600 kms" 1 . The first is seen in Siiv, Civ, 
N v, Si II, and C II. The second is seen in C iv, Si IV, and Si II. 
These two HV absorption-line complexes are shown in sep- 
arate columns in Figure 5. The Ovi data are too blended 
in the HV components to be of use. A full s tudy of the 
LIVES spectrum of GRB 021004 is presented by Fiore et al. 
(2005), and possible origins of the HV components (e.g. a 
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Table 7. GRB 071031 z=2.6922: Voigt profile fits 



Ion 




b 


log N 


log N (total) 
(N in cm" 2 ) 




(kms -1 ) 


(kms -1 ) 


(N in cm -2 ) 


Ovi 


-560±3 


18±3 


14.52±0.08 


14.77±0.08 




-505±2 


14±2 


14.41±0.06 




Ovi 


-118±7 


13±12 


13.79±0.42 


W15.14* 




-81±15 


* 


* 






-31±12 


* 


* 






6±14 


* 


* 






31±15 


* 


* 




Nv 


25±2 


11±2 


14.41±0.15 


14.41±0.15 


Civ 


-566±2 


* 


* 


?sl4.40* 




-559±5 


10±3 


13.06±0.27 






-515±3 


* 


* 






-505±4 


7±3 


13.67±0.24 






-368±2 


12±2 


13.93±0.03 






-325±2 


16±3 


13.14±0.05 




Civ 


-212±2 


4±3 


12.47±0.11 


wl5.ll* 




-147±2 


11±2 


13.13±0.06 






-111±2 


14±2 


13.96±0.08 






-70±6 


* 


* 






-35±14 


33±30 


14.43±0.44 






8±7 


22±5 


14.69±0.22 






39±5 


8±3 


14.58±0.36 




Si iv 


-367±2 


6±2 


12.33±0.08 


12.33±0.08 


Si iv 


-109±2 


15±2 


13.28±0.03 


15.93±0.51 




-63±2 


19±3 


13.61±0.07 






-40±2 


4±16 


15.89±0.48 






- 1±3 


29±3 


14.22±0.07 






29±3 


11±2 


14.65±0.24 




Cm 


-560±2 


15±3 


13.20±0.07 


«14.47* 




-507±2 


* 








-367±2 


* 







NOTE: the symbol * denotes a saturated component; no reli- 
able column density or line width can be extracted for the com- 
ponent, and the total column density in that ion is uncertain. 
NOTE: Total column densities (summed over components) are 
reported separately for the HV and LV absorption. NOTE: O VI 
is only detected at — 118 kms -1 in A1031 (A1037 is blended), so 
the detection in this component is uncertain. 

Table 8. GRB 080310 z=2.4274: Voigt profile fits 



Ion 


vo 

(kms -1 ) 


b 

(kms -1 ) 


log TV 
(N in cm -2 ) 


log JV (total) 
(N in cm -2 ) 


Ovi 


-137±4 


25±5 


14.18±0.08 


15.12±0.09 




-30±3 


42±4 


15.07±0.08 




Nv 


-49±2 


18±3 


13.75±0.05 


14.03±0.05 




1±2 


19±3 


13.71±0.05 




Civ 


-1414±2 


4±2 


12.89±0.10 


13.45±0.08 




-1388±4 


33±5 


13.31±0.07 




Civ 


-222±16 


22±14 


13.03±0.41 


^16.66* 




-187±4 


17±10 


13.36±0.31 






-160±3 


11±4 


13.33±0.25 






-138±2 


5±2 


13.44±0.36 






-66±50 










-10±50 


* 


* 






69±2 


5±2 


12.88±0.08 




Si iv 


-136±2 


4±2 


12.53±0.05 


«15.00* 




-55±50 




* 






-29±19 




* 






9±18 









NOTE: the symbol * denotes a saturated component; no reli- 
able column density or line width can be extracted for the com- 
ponent, and the total column density in that ion is uncertain. 
NOTE: Total column densities (summed over components) are 
reported separately for the HV and LV absorption. 



series of shel ls arou n d the progenito r star) a re discussed in 
iMoller et al.l (12001) . iMirabal et all (|2003l) . iSchaefer et all 
(<2003f ). and lLazzati et all (|2006f ). The velocity difference 
between the absorber at zgrb and the HV absorber near 
—500 kms -1 is equal to the velocity difference between the 
two lines of the C iv doublet, i.e. line-loc king is observed in 
Civ (first noticed bv lMoller et"aIll2002r Fl . 

3.2. GRB 050730 z GRB =3.9686 

High-ion absorption is seen in Ovi, Nv, Civ, Siiv, Siv, 
and Svi centered at zgrb- All of these lines except Nv 
are saturated. The Nv profile shows two components sep- 
arated by 27±5 kms -1 , each with fairly narrow b- values 
(7±3 and 14±8 kms -1 ); we fit the S vi absorption with one 
component with 6=17±2 kms -1 . Siv* A1072.973 is not de- 
tected at zqrb; we measure log V(S iv*)<13.70 (3er) over 
the velocity range —25 to 20 kms -1 , in which we measure 
logiV(Siv)=14.61±0.15, so that V(Siv*)/V(Siv)<0.12. 
Note that this limit applies only to the strongest com- 
ponent of absorption. The Civ, Siiv, Ovi, and Cn pro- 
files show a clear negative- velocity absorption-line wing ex- 
tending from —70 to —145 kms -1 . Finally, a HV absorber 
at —1550 kms -1 is seen with at least five Si iv and Civ 
comp onents spread over «200 kms -1 (see also IChen et al.l 
I2007D . This HV absorption-line complex is shown in Figure 
5. This GRB shows an extremely high Hi column of 
logiV(Hl)=22.15, and is the highest redshift GRB in our 
sample. A full study of the neutral and excited-state lines 
at zqrb in the UVES s pectrum of GRB 050730 is presented 
bv lD'Elia et all (|2007l) . 

3.3. GRB 050820 z GRB =2.6147 

Nv shows a single narrow component at —109 kms -1 rel- 
ative to zq RB with £>=9±3 kms -1 . This Nv-zgrb offset is 
the lar gest among our seven sp ectra, and has been noticed 
before ( Prochaska et al . 2008b), but the formal significance 



of the A1238.821 and A1242.804 detections is only 3.1cr and 
2.0cr respectively. The C iv and Siiv profiles exhibit at least 
ten components spread over 400 kms -1 in velocity width, 
with a Civ/Siiv ratio that changes dramatically between 
components. Civ A1550.781 is blended so we only display 
Civ A1548.204. A similar velocity structure is observed at 
lower S/N in Aim A1862.789. Svi A933.378 shows three 
components aligned with three components seen in C IV, 
strongly suggesting the genuine presence of S vi. However, 
Svi A944.523 is affected by low S/N and possible blend- 
ing, so we treat the S vi detection as marginal. The O vi 
profiles are highly blended and/or saturated, except for the 
range -140 and -200 kms" 1 in Ovi A1031.926, where the 
profile closely follows the shape of the C iv and Si IV pro- 
file. Therefore while we are not able to measure the O vi 
absorption in this case, we know that it is present. Siv 
A1062.662 shows a profile that broadly follows the Siiv. 
The Siv* A1072.973 profile is difficult to interpret. An ab- 
sorption feature at 15 kms -1 does not correspond to any 
component in any other high ion, so we identify it as a 
blend. However, we cannot rule out the presence of genuine 



° Line-locking is understood as occurring when a cloud is ac- 
celerated by radiation pressure in a given doublet line until the 
wavelength shifts int o the shadow of the ot her doublet line cast 
by a separate cloud (Vilkoviskii ct al. 1999). 
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Siv* absorption between —50 and kms -1 . Finally a two- 
component HV absorber at —3 660 kms -1 relative to zgrb 
is seen in Civ, but with no detectable Siiv, Cn, or Sin 
(see Figure 5). 

3.4. GRB 050922C z GRB =2.1990 

A single, strong component in Ovi, Nv, Civ and Siiv 
is detected at 65 kms -1 relative to zqrb- Aim is also de- 
tected in this component in both A1854.716 and A1862.789. 
Our best-fit to the O VI line has b=2A±2 kms -1 ; for the N v 
line, we find 6=16±3 kms -1 . C iv and Siiv show a complex 
series of at least six components spread over «280 kms -1 , 
but with no evidence for a negative- velocity wing. The S iv 
and S VI data are blended and not shown on Figure 3. This 
GRB shows a high logiV(Hi) value of 21.6 0, and a low 
metallicity of -2.0 (|Prochaska et all l2008bl ). This is the 
only GRB afterglow spectrum in our sample in which no 
HV C iv components are detected. 

3.5. GRB 060607 z GRB =3.0749 

The high ions in the GRB 060607 spectrum behave differ- 
ently than in the other afterglow spectra. The Civ, Siiv, 
and Si in profiles show a four-component absorber centered 
at zgrb and covering only wlOO kms -1 of total width. 
However, Nv is not detected at ^grb; we derive a 3a up- 
per limit of log A^(N v)<12.70 in the velocity range —140 
to 50 kms -1 , where Civ absorption is seen. The O vi, Siv, 
and Svi data are blended near zqrb and are of no use. 
Aim shows a non-detection. The Civ profiles also show 
a weak (low optical depth) wing extending from —60 to 
— 120 kms -1 . This GRB has an extremely low H i column 
density, log iV(Hl)=16.8 (jProchaska et al.l l2008bh . 3.5 or- 
ders of magnitude lower than the DLA columns seen in five 
of the other six cases. A HV absorber seen in C iv and Siiv 
is centered at —1850 kms -1 relative to zqrb, with nearby 
low- ion absorption present in Sin and Cn (Figure 5). 

We note that a sub-DLA at z=2.9372 is detected in the 
GRB 0606 07 spectrum with log jV(Hl)=19.50±0.05 (first 
noticed bv lLedoux et al.ll2006f ). This sub-DLA is separated 
from the nominal GRB redshift of 3.0749 by Az=0.1377, or 
10 300 kms -1 , and shows strong, multi-component absorp- 
tion in O vi, C iv, and Si iv. We considered the possibility 
that the host galaxy of GRB 060607 is actually located at 
this redshift (z=2.9372), which would imply its ISM prop- 
erties are more similar to those of the host galaxy ISM of 
the other bursts in our sample. However, this hypothesis 
is unlikely to be correct since (a) Lyman-a forest lines are 
detected up to the redshift of 3.0749, suggesting that zqrb 
really is this high, and (b) this scenario would require the 
absorber observed at z=3.0749 to be in front of the host 
galaxy yet moving toward it at over 10 000 kms -1 , which 
is not easy to explain dynamically. We conclude that the 
sub-DLA at 2.9372 traces a foreground galaxy unrelated to 
the GRB. 

3.6. GRB 071031 z GRB =2.6922 

A strong component in Ovi, Nv, Civ, and Siiv is de- 
tected at zqrb- The least-saturated high-ion line, Nv, 
shows a narrow line width of 6(Nv) = ll±2 kms -1 , imply- 
ing logT < 5.08 for the gas containing the N +4 ions, re- 



gardless of ionization mechanism. No information is avail- 
able on the line widths of Civ and Siiv at zq RB be- 
cause the absorption in these lines is saturated. Whereas 
Svi A933.378 is blended, Svi A944.523 appears to show 
absorption at zqrb, with a profile that matches that of 
Ovi. However, because of the low S/N of the data at 
these wavelengths, we do not attempt to derive a Svi 
column density. We measure log N (S IV*)<13.87 (3a) over 
the velocity range to 50 kms -1 , but complete satura- 
tion (and possible blending) prevents a reliable measure- 
ment of logiV(Srv) in this velocity range. Absorption in 
several fine-structure lines of Fe, Ni, O, Si, and C is de- 
tected at zqrb in the GRB 071031 afterglow spectrum: in 
Figure 4 we show a line (Fen***** A1559.085) from the 
fifth excited level ( 4 F g / 2 ) of Fen. We derive an Hi column 
density log iV(Hi)=22.15±0.05 and a metallicity [Zn/H] of 
— 1.72±0.06, where zinc is chosen since it is typically unde- 
pleted onto interstellar dust grains. 

A notable feature of the C iv, Si iv, S iv, and (at lower 
S/N) the Ovi profiles is the wing extending blueward 
from kms -1 down to —160 kms -1 . A component at 
— 115 kms -1 is superimposed on the wing. The Civ and 
Siiv profiles follow one another closely, suggesting these 
two ions form in the same physical regions of gas. However, 
the slope dN a (v)/dv is shallower for Ovi than for Civ and 
Siiv in the range —90 to —50 kms -1 (see §4.3). 

Three strong components of HV high-ion absorption 
are seen (Figure 5), at velocities of —560, —510, and —370 
kms -1 . The —560 and —510 kms -1 components are seen in 
Civ A1548.204, Ovi A1031.926, Ovi A1037.617, and Cm 
A977.020, but are blended in Civ A1550.781 and absent in 
Si iv. The -370 km s -1 component is detected in C in, C iv, 
and Siiv, but the Ovi profiles are fully saturated at this 
velocity. The N v data show no significant detection in the 
HV gas, though the S/N near the Nv lines is too low to 
place a strong constraint on 7V(N v). Neutral-phase absorp- 
tion is not seen in the three HV components, even in the 
sensitive (high /-value) lines O i A1302.169, Sin A1260.422, 
and Cn A1334.532, with the single exception of a detection 
of weak Cn at -370 kms -1 . The V(Civ)/A(Siiv) ratios 
change dramatically between the various components, from 
4.8±1.0 at -110 kms -1 , to 40±9 at -370 kms" 1 , >140 
(3a) at -510 kms -1 and >110 (3a) at -560 kms -1 . 

3.6.1. Time-series observations of GRB 071031 

The combined UVES spectrum of GRB 071031 (shown 
in Figure 4) was formed by co-adding four useful indi- 
vidual exposures, with exposure times of 300, 600, 1200, 
and 2400 s (the sequence of increasing exposure times was 
chosen to compensate for the fading of the afterglow). 
Analysis of these individual frames allows us to look for 
time-variation in the high-ion absorption. We look for vari- 
ation in the Civ and Siiv profiles, since these lines fall 
in parts of the spectrum where the S/N is highest (a Nv 
time-series is unavailable). We focus on the HV compo- 
nents, since the high-ion absorption at zgrb is saturated, 
making time- variation at that redshift difficult to detect. 
The time-series spectra are shown in Figure 6. No evi- 
dence for any evolution in the C IV or Si IV profiles is seen. 
The lack of variation can be quantified by measuring the 
equivalent width at each observing epoch. The total equiv- 
alent width and la error measurements in the range —600 
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GRB071031 C IV 1548 



GRB071031 Si IV 1393 




Fig. 6. Time-series observations of C iv and Si IV absorp- 
tion in the GRB 071031 afterglow spectra. Each row shows 
a different observing epoch, with the UT of the start of 
each exposure annotated with the exposure duration on 
the right panel. These data have been binned by four pix- 
els to improve the S/N. No evidence for time- variation in 
the high-ion profiles is seen. 



to —330 kms 1 (over three HV components) are 266±88, 
327±49, 233±36, and 306±30 mA in the four exposures re- 
spectively, consistent with no variation in the absorption 
with time. Thus the behavior of the high-ion lines in these 
HV absorbers is in contrast to the time-variation in low- 
ion and fine-structure absorp t ion li nes at zqrb reported 
by Dessauges-Zavadskv et~ahl (|2006l ) and IVreeswiik et al.l 
(200^ _ 



3.7. GRB 080310 z GRB =2.4274 

Strong high-ion absorption is seen in O VI, N v, C iv, S IV, 
and Si iv centered within 50 kms -1 of zgrb- The Nv pro- 
file shows two components separated by 50±3 kms™ 1 and 
with 6-values of 18±3 and 19±3 kms™ 1 . The Svi profiles 
are suggestive of absorption at zgrb, but with very low S/N 
and an inconsistency between the absorption in A933.378 
and in A944.523. For this reason, we do not attempt to 
derive a Svi column density. Similarly the Siv A1062.662 
profile is too noisy to measure a reliable column, but we 
can measure a lower limit of log N(S iv)>14.77 over the ve- 
locity range —70 to 30 kms™ 1 . Combining this with the 
constraint logiV(S iv*)<14.38 (3er) from the non-detection 
of Siv* A1072.973, we find JV(Srv*)/JV(Srv)<0.41. The 
Civ, Si iv, and Ovi profiles show a clear negative- velocity 
absorption-line wing extending from —120 to —260 kms™ 1 . 
This is the most extended absorption-line wing in our sam- 
ple. Finally, a HV absorber at —1400 kms™ 1 is seen (see 
Figure 5) with two C iv components, one of which is very 
narrow (6=4±2 kms™ 1 ). For this GRB-DLA we derive an 
Hi column density log A r (Hi)=18.80±0.10 and a metallicity 
[O/H] of — 1.39±0.10. No ionization correction was applied 
since a charge- exchange r eaction closely links O I and Hi, 
implying [O i/H l]«[0/H] dField k Steigmadll97l . 



4. Discussion 

The six high-ion absorption lines reported in this paper 
at velocities near zgrb are Si IV, Siv, Civ, Svi, Nv, and 
Ovi, which trace the Si+ 3 , S+ 3 , C+ 3 , S+ 5 , N+ 4 , and 0+ 5 
ions, requiring energies of 33.5, 34.8, 47.9, 72.6, 77.5, and 



113.9 eV for their creation, respectively (|Mortonll2003[ ). A 
galactic ionizing spectrum, dominated by the integrated ra- 
diation from O and B s tars, drops strongly above 54 eV, 
the Hell ionization edge (|Bregman fc Harrin gton 1986ft . So 
whereas the ions Si +3 , S +3 , and C +3 can be photoionized in 
the ISM by starlight, the ions S+ 5 , N+ 4 , and 0+ 5 cannot. 
We expect this to remain true in GRB host galaxies, even 
though they have be en shown to harbor s ignificant numbers 
of Wolf-Rayet stars (|Hammer et al.|[2006h . since Wolf-Raye t 
spectra also show a strong break at 54 eV ((Crowthei][2007). 
Instead, the detection of the three ions S +5 , N+ 4 , and +5 
implies the presence of either a hard, non-stellar radiation 
source or hot, collisionally ionized gas. 

In the seven GRB afterglow spectra in our sample, the 
median and standard deviation of logiV(Ovi) at zqrb is 
15.12±0.54. For comparison, the median logA(Ovi) mea- 
sured in 100 sight lines passing through the Milky Way 
halo is 14.38 (|Wakker et alj|2003t ISavage et al.ll2003l ). ris- 
ing to 14.80 w hen including the con tribution from high- 
velocity clouds (jSembach et al.|[2~003f h Turning to Nv, we 
report a median log iV(N v)=14.03±0.49 in the seven GRB 
spectra, whereas t he median log N (N v) in 32 Galactic halo 
sight lines is 13.45 (jlndebetouw k Shulfeool . High-ion ab- 
sorption in Nv and Ovi has also been observed in z=2-3 
PL A g alaxies, with a median log N(0 vi) = 14.77 (|Fox et al.l 
12007a) ). and 76 per cent of DLAs showing log iV(N v)<13.50 
(Fox et al. 2008, in preparation). Thus the column densi- 
ties of Ovi and Nv seen at Zgrb a re , m the mean, higher 
than those seen in the halo of the Milky Way and those 
in DLAs. The highly ionized Galactic ISM is a complex 
environment, with many physical processes contributing to 
the production of the high ions, including conductive inter- 
faces, turbulent mixing layers, shocks, e xpanding supernova 
remnants, and galactic fountain flows (ISavage et all [2 003; 
I Zsargo et all 120031: Ihidebetouw fc Shulll 12004 iBowenet all 
2008), and there is no reason why the ISM in the host 
galaxies of GRBs should be any less complex. 

With this in mind, we now discuss the various categories 
of high- ion absorption seen near zgrb- Throughout this sec- 
tion, we refer the reader to Table 9, which summarizes the 
high-ion to high-ion column density ratios measured in the 
GRB afterglow spectra, and compares them to the ratios 
measured in various other astrophysical environments. 

4.1. Strong absorption at zgrb 

We always detect strong high-ion absorption exactly at 
zgrbEL in the form of saturated Ovi, Civ, and Siiv com- 
ponents, and (where the data are unblended) also in SlV 
and Svi. Nv is present in 6/7 cases, with the advantage of 
being less saturated than the other high-ion lines. Because 
of this, Nv is the best line to investigate the line width 
and optical depth in the strong absorption components. In 
the six cases where N V is detected, we fit eight N v compo- 
nents, with b- values of 7±3, 9±3, 11±2, 14±8, 16±3, 18±3, 
19±3 and 28±3 kms™ 1 . The median of these eight values 



In two cases the velocity centroid of the strongest high-ion 
absorption component differs from the nominal GRB redshift 
(measured from the low-ion lines) by more than 10 kms™ 1 : for 
GRB 050922C the offset is 65 kms™ 1 , and for GRB 071031 the 
offset is 25 kms™ 1 . However, since the GRB redshifts are not 
known to better than a precision of several tens of km s~ 1 , these 
offsets are not significant, and we treat the strong high-ion ab- 
sorption as arising exactly at zgrb. 
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Table 9. High-ion column density ratios 



GRB^ 


W(Oiv) 

N(0 vi) 


JV(Oiv) 
JV(Silv) 


N (N v) 
Af(0 vi) 


At z(GRB) 








050922C 






0.21±0.04 


060607 




7.6±1.2 




Blueshifted Wings 








050730 


0.40±0.25 


3.7±0.9 


<0.09 


UoUoU ( 




>lo 




U ( ±UO± 




9-t-n q 
o.zztu.y 


<.U.ll 


080^1 


u.uuznu. J. 1 


25±4 




I I V/ pnninnnpnts 

XX V LUU1 UUllCllUo 








021004 -2900 kms" 1 


* 


5.1±3.1 


* 


050730 —1550 kms -1 


* 


3.7±0.4 


* 


050820 —3660 kms -1 


* 


>24 


* 


060607 — 1850 kms -1 


* 


3.2±0.5 


* 


071031 -560 kms" 1 


0.19±0.06 


>110 


<0.06 


071031 -510 kms -1 


0.35±0.10 


>140 


<0.07 


071031 -370 kms" 1 




40±9 


* 


080310 -1400 kms -1 


* 


>36 




Galactic Halo 6 


0.60±0.47 


3.46±1.09 


0.12±0.07 


LMC_ 


<0.09-0.59 


1.5-2.5 


<0.33 


DLA~ z=2-3^ 


2.1±0.7 


5.4±2.2 




IGM, 2=2-4 ~ 


0.35 c 


16±4 / 


<0.23 9 



a We report the observed column density ratios in different 
velocity regions: in the strong components seen at zgrb, in the 
blueshifted wings, and in the HV absorbers. All limits are 3a. In 
cases marked *, saturation, blending, or low S/N prevents the 
ratio from being measured reliably. 

b M ean±lq in 16 Milky Way halo sight lines (|Zsarg6 et al] 
120031 ). 

c Range obser ved in LMC gas (not including H II regions) in 
four sight lines ()Lehner fc Howld 120071). 

d Meanila in 10 DLAs (|Fox et alj|2007aP ). 

e Median value among four IGM absorbers measured by 
ISimcoe et all (|2006l ). 

f Medianilcr in 188 IGM absorbers (jBoksenberg et al.ll2003l). 

9 M easured in composite QSO spectrum by ILu fc Savage! 
lfl993T) . 



is 16 kms -1 . A Nv component formed in gas at 200 000 K, 
the temperature at which the production o f Nv by colli- 
sions is maximized ()Gnat fc Sternbe rg 2007), has a b- value 
of 15.4 kms -1 . Therefore, we are unable to rule out colli- 
sional ionization for the gas traced by Nv. 



4.1.1. Photoionized circumburst gas? 

Noticing that the Nv absorption at zqrb tends to be 
stronger and narrower (in total velocity width) than the 
Nv a bsorption observed in th e Milky Way ISM and in 
DLAs, iProchaska et al] (|2008bl) have recently argued that 
the Nv absorption components at zgrb do not trace ISM 
material in the halo of the host galaxy, but rather arise in 
circumburst gas in the immediate vicinity of the GRB. Our 
high-resolution UVES dataset, covering a wider range of 
ionization states than has been observed before, allows us 
to investigate this intriguing idea. 

Our data are qualitatively consistent with the circum- 
burst hypothesis in the cases of GRBs 021004, 050922C, 
and 071031, where we see strong, single components of Nv 
at zgrb, each with log N(N v)>14.0, aligned with other 
high ions including Ovi with log A(0 vi)>14.5. For these 



cases, our new measurements of the O VI, S IV, and S vi 
column densities at zgrb (listed in Tables 2, 4, and 7) 
could be used to constrain photoionization models of the 
circumburst region. However, in the remaining four spectra 
in our sample (GRBs 050730, 050820, 060607, and 080310) 
the situation is more complex. In each of these four cases, 
logA(Nv)<14.1, and in two cases log N(N v)<13.2, i.e the 
observed N v is not in excess of what is expected from the 
host galaxy's ISM. In two of these cases (GRBs 050730 and 
080310) the N v absorption profiles show two components 
separated by a few tens of kms -1 - this is not the signa- 
ture of a single burst. In another case (GRB 050820) there 
is an offset of 110 kms -1 between a single Nv component 
and the strongest absorption in the other high ions. And 
finally, in the case of GRB 060607, no Nv is detected at 
all. In addition to these case-by-case details, we repeat that 
the b- values of the detected N V absorption components at 
zgrb are not particularly narrow, so photoionization is not 
required by the line widths (though it is not ruled out ei- 
ther). 

4.1.2. Total ionized column density 

We can further investigate the nature of the strong com- 
ponents seen at zqrb by calculating the total column den- 
sity of ionized hydrogen JV(hot Hn)Q contained in these 
components. To do this we have to correct for ioniza- 
tion and metallicity, so the calculation is only possible 
for cases where the nitrogen abundance has been mea- 
sured in the DLA at zgrb- This is true for three GRBs 
in our sample: [N/H]=-3.16±0.10 for GRB 050730, > 
-1.35 for GRB 050820 , and < -4.09 for GRB 050922C 
()Prochaska et alj|2007bl ). The relationship between the Nv 
and hot H II column densities can be written in the follow- 
ing way: 



A (hot Hil) 



JV(Nv) (N/H) n 
(Nv/N)(N/H) n (N/H)i ' 



(1) 



where (N/H)i denotes the nitrogen abundance in the ion- 
ized gas and (N/H) n denotes the nitrogen abundance in the 
neutral gas (which is measured). We take the solar nitro- 
gen abundance of (N/H) = 10~ 422 from iGrevesse et al] 
d2007h . and then (N/H) n = 1O[ N / H 1(N/H) . If the gas 
is collisionally ionized, as would be appropriate for Nv 
in the host galaxy ISM, Nv/N<0.25 at all temperatures 
(jGnat fc SternberglteOOTli 5 ! . Therefore we adopt 0.25 as the 
maximum allowed ionization fraction, corresponding to a 
minimum ionization correction of a factor of 4. If the ion- 
ized gas has the same metallicity as the neutral gas, so that 
the term (N/H) n /(N/H)i in Equation 1 is equal to one, then 
AT (hot Hn) is: 

>21.8 for GRB 050730 [which has log A(Hi)=22.15], 
>19.3 for GRB 050820 [which has log A(Hi)=21.00], and 
>22.7 for GRB 050922C [which has log A(Hi)=21.55], 
where the results are lower limits since (a) we have used 
the smallest allowed ionization correction, and (b) they 
only apply to the Nv-bearing gas, and Hil may exist 



7 Here we use the word "hot" simply to label the N v-bearing 
gas; the upper limit on the temperature is 215,000 K for the 
median Nv b- value of 16 kms -1 . 

8 This is true f or either equilibrium or n on-equilibrium models. 
However, in the [Prochaska et al. (2008b) GRB photoionization 
models, a Nv ionization fraction as high as 0.6 is predicted. 
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Fig. 7. Analysis of the S IV level populations at zqrb to- 
ward GRB 050730. The solid and dashed lines show the 
predictions from UV photo-excitation models for the time- 
variation of the Siv and Siv* column densities (see text 
for details). The data points show the observed 7V(Srv) 
and upper limit to AT(Siv*) in the UVES spectrum, plot- 
ted at the x-position corresponding to the mid-observation 
epoch (five hours after the trigger). The solid lines show the 
model with the best-fit distance of 418±9 pc. The dashed 
lines show a model with the distance d fixed at 200 pc, for 
illustration. In this case N(S IV*) would be 0.5 dex stronger 
than observed, so this model is clearly ruled out. The result 
that the S iv-absorbing gas must lie at d > 400 pc implies 
an interstellar rather than circumburst origin. 



at other temperatures as well. Comparing the values of 
7V(hot Hn) and A^(Hi), we find the hot-ionized-to- neutral 
ratio iV(hot Hll)/iV(Hl) takes values of >0.4 and >0.02 
for GRBs 050730 and 050820, and >13 for GRB 050922C. 
Whereas the first two values are modest, the GRB 050922C 
value is extremely (and implausibly) high - it requires over 
ten times as much mass in hot ionized gas in a single com- 
ponent as there is in neutral gas in the entire host galaxy. 
However, if the nitrogen abundance in the ionized gas was 
higher than in the neutral gas, so that (N/H) n /(N/H)j < 1, 
then solutions with lower iV(hot Hn) would be possible. 
Thus our calculation of the total ionized column density 
at zgrb suggests that the metallicity in the highly ionized 
gas is higher than the metallicity in the neutral gas (par- 
ticularly for GRB 050922C); however, this calculation is 
inconclusive in determining whether the gas is circumburst 
or interstellar. 

4.1.3. Photo-excitation modeling of iV(Siv) and N(Sw*) 

Another method for distinguishing between circumburst 
gas and more distant interstellar gas is to look for ab- 
sorption from fine-structure levels. In circumburst gas, UV 



pumping by GRB photons will populate the excited elec- 
tronic levels, which decay and cascade to populate the 
fine-structure levels of the ground state (IProchaska et al.l 
120061: IChen et al.ll2007t IVreeswiik et~"aT]|2007l ). Because of 
the rapid fading of the afterglow, the level populations 
evolve in a non-equilibrium, time-dependent manner. Our 
data allow us to investigate this process in the S IV ion, 
since in four GRB spectra the combination of a detection 
of Siv A1062.662 and a non-detection of Siv* A1072.973 
allows us to place an upper limit on the SlV*/SlV ratio 
at zgrb- The Einstein A- value of the forbidden transi- 
tion between the ground state of S IV and its fine-structure 
level is 7.7 x 10" 3 s _1 (taken from NIST: iRalchenko et al.l 
^ In other words, in the absence of stimulated photo- 



and collisional de-excitation, the fine-structure level probed 
by the Siv* 1072.973 transition will decay in 1/A=130s. 
In fully ionized gas (as expected near GRBs) there is no 
Lyman-a; opacity and so the S iv lines (which lie shortward 
of Lyman-a) will not be shielded from UV radiation. For 
GRBs 021004, 050730, and 080310 we measure 3cr limits 
to the N(Sw*)/N(Sw) ratio of <0.08, <0.12, and <0.41 
respectively (for the other four afterglows we are unable to 
constrain this ratio due to blending). 

Fo llowing the technique described in IVreeswiik et alj 
(l2007l) . we can model the Siv excitation level in gas 
subject to UV radiation from the GRB, assuming that 
photo-excitation is the dominant excitation process. We 
chose to model the S IV level population measured toward 
GRB 050730, because this case has the strongest constraint 
on N(Siv*). The model has five free parameters: the dis- 
tance d between the GRB and the absorbing gas (assuming 
all the S iv arises at the same distance) , the total S iv col- 
umn density N to t (equivalent to the pre-burst ground-state 
S iv column density) , the afterglow spectral slope (3 (where 
F v oc v^), the Doppler 6-parameter of the absorbing gas, 
and the rest-frame time when the calculation begins to (i.e. 
the time when the afterglow photons begin to be tracked in 
the code). For a given set of these parameters, the model 
predicts the S iv and S iv* column densities as a function 
of time (both in the observer frame and in the rest frame). 
Since the time of mid-observation of GRB 050730 is known, 
five hours (observer frame) after the trigger, predicted val- 
ues for N(Sw) and A^Siv*) can be extracted from each 
model. 

The parameters in our GRB 050730 model were 
determined as follows. N to t was chosen to equal 
7V(Siv)+7V(Siv*), assuming the actual iV(Siv*) is equal 
to the measured upper limit (if the actual value is lower, 
our conclusions are strengthened; see below) . The value for 
(3 (—0.56) was determined using the observed light curve 
for this afterglow. The parameter to is poorly constrained, 
but our runs showed that the S IV level populations are 
very insensitive to to for values between and 300 s; we 
used a value io=10s (see IVreeswiik et al"1l2007l ). Finally a 
6- value of 10 kms -1 was chosen as a typical value for in- 
terstellar gas, and indeed matches (within the l-a error) 
the measured b- value for the Siv component at kms -1 
toward GRB 050730. This leaves the distance d as the one 
genuinely floating parameter in the model; a grid of mod- 
els at different distance d was generated, and the predicted 
and observed N(S iv) and N(S iv*) were compared using a 
chi-squared minimization routine. 



See http://physics.nist.gov/PhysRefData/ASD/index.html 
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The model results are shown in Figure 7. The best-fit 
distance d is 418±9 pc. Because we assumed that the ac- 
tual S iv* column density was equal to the measured upper 
limit, this represents a strong lower limit of 400 pc (2a) on 
the distance of the absorbing gas from the GRB. In other 
words, if the gas was at d < 400 pc, we would have seen 
Siv* absorption, whereas we clearly do not. This repre- 
sents an important result, since it implies that the Siv ab- 
sorption has an interstellar rather than circumburst origin. 
The distance limit on the S iv absorption does not necessar- 
ily have to apply to the Nv absorption (or the other high 
ions). However, the Nv and Siv profiles at zqrb toward 
GRB 050730 are similar, both showing the component near 
krns" 1 and a weaker component at 25 kms -1 , suggest- 
ing that the two ions are co-spatial and hence that the Nv 
absorption also arises at >400 pc from the GRB. 



4.1.4. Origin of strong absorption 

In summary, we find that while the circumburst hypoth- 
esis explains the observation that in three afterglow spec- 
tra, strong single-component N v is aligned with saturated 
absorption in the other high ions, it is challenged by the 
Nv profiles in the other four datasets, it is not strictly re- 
quired by the Nv line widths, and the expected high-ion 
fine-structure lines and time-variable column densities are 
yet to be observed. Furthermore, in one case (GRB 050730) 
we are able to place a lower limit on the distance to the 
high-ion absorbing gas of 400 pc, using the simultaneous 
detection of Siv and non-detection of Siv*. Therefore, 
based on our current dataset, we cannot rule out an inter- 
stellar explanation for the strong high-ion components at 
z grb, noting that since the higher iV(Hi) values in GRB- 
DLAs than in QSO-DLAs suggest that GRB sightlines pass 
through the inner star-forming regions of t he host galaxies 
(|Prochaska et alj|2008at iFvnbo et al.ll2008h . one would ex- 
pect high column densities of interstellar plasma at zqrb , 
so the strength of the Nv absorption does not, in our view, 
argue against an interstellar origin. We also note that a 
strong, apparently narrow (though unresolved) Nv com- 
ponent is observed in t he spectrum of LB G cB58 at the 
galaxy redshift (z—2.73\ iPettini et al.ll2002t ). showing that 
such components can be formed even in the absence of a 
GRB. Time-series observations of Nv in GRB afterglow 
spectra, which unfortunately are not available in the cur- 
rent dataset, would help to resolve this issue. A significant 
detection of time- va riation in the Ny colum n density at 
zgrb (see models bv lProchaska et al"1l200 8b') would repre- 
sent the discovery of gas that is unambiguously close to the 
GRB. 



4.2. Multi-component absorption in Civ and Siiv 

The multi-component structure seen in C iv and Si iv 
spreading over several hundred kms -1 around zqrb is 
not seen in the other high ions. However, the C iv and 
Si iv profiles lie out of the Lyman-a forest, so are free 
from blending and show higher S/N than Ovi or Nv, 
and it is unclear how the Ovi and Nv profiles would ap- 
pear if observed at the same S/N. These complex, multi- 
component Si iv and C iv profiles are reminiscent of the 
profiles of these i ons in damped Lyma n -q (DLA) absorbers 
(|Lu et alJ 119961: iLedoux et all 119981 : IWolfe fc Prochaskal 



12001 IFox et al] 120071* iLehner et al.ll200l . which repre- 
sent galaxy halos seen in absorption toward a background 
source. The median (mean) value of the total Civ veloc- 
ity width_Jn 74T2LAs and sub-DLAs at z=2-3 reported 
bv IFox et all (|2007bD is 255 kms -1 (342 kms" 1 ), where 
the width is the total observed range of absorption, re- 
gardless of optical depth. In the seven GRB spectra in 
our current sample, the median (mean) total C IV veloc- 
ity width measured in an identical manner is 280 kms -1 
(320 kms -1 ), similar to the DLA value (note we have ex- 
cluded the HV components, which we treat separately). A 
precise comparison of the median C iv column density be- 
tween the GRB sample and the DLA sample is not possible, 
because these measurements are affected by saturation, but 
we can state that log iV (Civ) at zgrb is ^15 in all seven 
GRB spectra, which is significantly stronger than the me- 
di an DLA value log jV(DLA C iv)=14.27±0.57 measured 
bv IFox et al] (|2007bl ). Nonetheless, the similarities in total 
C IV line width between GRB and DLAs support a galactic 
origin for the multiple components seen in C IV and Si iv in 
a velocity range of several hundred kms -1 around zgrb- 



4.3. Negative-velocity (blueshifted) absorption- line wings 

One notable feature present in four of the seven GRB 
afterglow spectra in our sample is a negative-velocity 
absorption- line wing seen in Civ, Siiv, and (at lower S/N) 
Ovi. These wings are characterized by relatively smooth, 
asymmetric absorption extending for 65-140 kms" 1 to neg- 
ative velocities relative to the strong absorption at zgrb- 
To explore the high-ion velocity structure in the wings, we 
present in Figure 8 apparent column density profiles of the 
high ions over the velocity of the wings, calculated using 
N a (v) — 3.768 x 10 14 T a (w)//A, where the apparent optical 
depth T a (v) = \n[F c (v)/F(v)], F(v) and F c (v) are the ob- 
served flux level and estimated continu um level as a func- 
tion of velocity, and where A is in A (Savage & Sc mbach] 
1 19911) . We follow the error treatment of lSembach fc Savage! 
(1992f). This form of display allows the velocity structure 
seen in each ion to be closely compared. In two cases we 
add the profile of Cn A1334.532, which also shows wing 
absorption. The high-ion column densities in the wings are 
summarized in Table 10. 

We find that the ClV and Siiv profiles track one an- 
other closely in three of the four wings, indicating that the 
two ions form in the same regions of gas. The exception is 
the very weak wing seen toward GRB 060607, in which Si iv 
is not detected, but Si in is seen. For the wings seen toward 
GRBs 050730 and 080310, the Ovi and Civ profiles share 
a common gradient dN a /dv. However, for GRB 071031, 
the Ovi absorption shows a shallower slope. We measure 
a 2V(Crv)/jV(Sirv) ratio of «3 in two of the four wings 
(GRBs 050730 and 071031), similar to the Galac t ic hal o 
average of 3.46±1.09 measured by IZsargo et all ((2003). 
This supports an interstellar origin for the absorption-line 
wings. In the wings toward GRBs 060607 and 080310 the 
iV(Clv)/JV(Silv) ratio is much higher, taking values of >13 
(3<t) and 25±4 respectively. 

O vi absorption-line wings with similar velocity extent 
have been observed at zero redshift in the ultraviolet spec- 
tra of quasars taken with the Far- Ultraviolet Spectr o scopic 



Explorer (FUSE) sat ellite (ISembach et al] l200ll [2 003: 



IFox et all 120051 l2006t ISavage et al] l2005at iKeenev et al I 
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Table 10. Properties of absorption in the blueshifted wings 



GRB 


ZGRB 


^min i ^ma>^ 


5t£ 


logiV(Ovi) 


logTV(Civ) 


logiV(Siiv) 


logiV(Nv) 






(kms" 1 ) 


(kmsT 1 ) 


(N in cm" 2 ) 


(N in cm" 2 ) 


(N in cm" 2 ) 


(N in cm" 2 ) 


050730 


3.9686 


-145,- 80 


65 


14.47±0.35 


13.79±0.03 


13.21±0.02 


<13.56 


060607 


3.0749 


-140,- 60 


80 




13.00±0.02 


<12.27 


<12.91 


071031 


2.6922 


-150,- 40 


110 


14.25±0.06 


14.45±0.09 


13.95±0.06 


<13.59 


080310 


2.4274 


-260,-120 


140 


14.31±0.08 


14.05±0.02 


12.65±0.06 


<13.27 



a On velocity scale relative to zgrb- 

b Velocity extent of wing <5u = « max -Umin, shaded light grey on Figures 2, 3, and 4. 




-140 -130 -120 -110 -100 -90 -80 -140 -120 -100 -80 -60 -40 




Fig. 8. High-ion apparent column density profiles in the four negative- velocity wings in our sample. Errors reflect both 
statistical noise and continuum placement uncertainties. The C iv and Si iv profiles have each been scaled by the factors 
annotated on the plot to allow for comparison of their shape with O vi. Only lower limits (arrows) can be presented for 
saturated pixels. The Ovi data have been rebinned by four pixels in the case of GRB 071031, eight pixels in the case of 
GRB 050922C, and five pixels in the other cases. Nv, SlV, and Svi are not detected in the wings. 



20060 iFox et alJ (|2006h measured eleven Galactic O vi 
wings, and report a mean and standard deviation O vi col- 
umn density of log7V(0 Vl)=13.85±0.45, with each wing ex- 
tending over ~100 kms" 1 ; these authors successfully mod- 
eled the O vi absorption-line wings as tracing Galactic out- 
flows moving under ballistic trajectories. Furthermore, a 
strong O vi wing covering a similar velocity range is seen 
in a sub -DLA (an absorb er with 19.0<logiV(H i)<20.3) at 
z=2.67 (|Fox et al.ll2007cj) . The Civ profile observed in the 
LBG cB58 at z=2.7 shows a clear absorption wing extcnd- 



1 Analogous C iv and Si IV wings are yet to be observed at 
zero redshift, but these two ions fall in the near-UV, requiring 
HST observations rather than FUSE, and the number of high- 
resolution QSO spectra taken with HST is small compared to 
the number of extragalactic sight-lines observed with FUSE. 



ing from —250 to —750 k ms" 1 (sa nine resol ution elements) 
in the galaxy rest-frame (|Pettini et al.ll2002T ) , and an asym- 
metric blueshifted absorption feature in Mg n has been re- 
po rted in the composit e spectrum of 1400 z w 1.4 galaxies 
bv lWeiner et alJ (|2008h . who also interpret it as a galactic 
outflow signature. 

Whereas the z=0 O VI wings are always detected at pos- 
itive velocities (redshifted) relative to the Local Standard 
of Rest (LSR) , the GRB wings are seen at negative veloci- 
ties (blueshifted) relative to 2grb • This behavior is strongly 
suggestive of galactic outflows: in the Milky Way case (at 
least for high-latitude and anti-center sight lines, where 
Galactic rotation effects are insignificant), outflowing gas 
appears redshifted relative to the LSR, whereas in the GRB 
case, outflowing gas moves toward us along the line-of-sight 
and hence appears blueshifted relative to the host galaxy. 
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Fig. 9. Montage of apparent column density profiles of O vi absorption-line wings measured in various environments: 
near zqrb in three GRB after glow spectra (th is paper, three left-most lines), at zero redshift in the Milky Way halo i n 
the UV spectra of three AGN (|Fox et al.ll2006l three right-most lines), and in a sub-DLA at z=2.6562 (|Fox et ai.ll2007cD . 
The Milky Way and sub-DLA cases have been inverted and velocity offsets have been added for ease of comparison, but 
the gradient of each profile is preserved. Linear fits are shown with color-coded dashed lines, and the gradient of each fit 
in units of 10 9 ions cm" 2 (kms -1 )" 2 is annotated with its lex error on the panel. In these diverse galactic environments 
sampling over twelve billion years, the Ovi absorption-line profiles can be used to trace and quantify galactic outflows. 



The possibility remains that some wing absorbers repre- 
sent chance alignments of unrelated components separated 
by small velocities, which are unresolved by the spectro- 
graph. However, we find this explanation unlikely for the 
strongest wings in our sample, seen toward GRBs 071031 
and 080310, since in these cases the column density in the 
wing increases monotonically with velocit\F T l in an interval 
of over 100 kms" 1 , and there is no reason why a series of 
random components should be aligned in this way. 

To illustrate the similarity in the wing absorption fea- 
tures observed at both low and high redshift, we show in 
Figure 9 a comparison of the O vi apparent column den- 
sity profiles in the zero redshift and GRB sight line wings. 
We have inverted the positive-velocity wings to compare 
their gradient with the wings observed in the GRB spec- 
tra. While the similarity in the overall shapes of the wing 
profiles is striking, the high-z wings in the GRB host galax- 
ies show higher O vi column densities (by a factor of «5) 
and steeper gradients dN a /dv (by a factor of 5-10) than 
the Milky Way wings, implying higher mass flow rates. 

4.4. High-velocity components at 500-5000 /cms" 1 

In six of our seven GRB spectra (all GRBs except 050922C), 
high-ion components at 500-5000 kms -1 relative to 
zgrb are observed (see Figure 5). Such HV components 
have been observed be fore in GRB 0208 1 3 dBarth et al. 
20031 ) GRB 021004 (iMoller et alj l2002t iMrabal et al. 



2003 : ISchaefer et alj [20031: iFiore et all 120051 : Eazzati et al 
20061) . GRB 030226 (iKlose et al.1 120017 GRB 050505 



(iBereer et al J 12006ft. and GRB 050730 (|D'Elia et all 120071: 
IChen et al.ll2007t ). IChen et all (|2007D report that C IV corn- 



In the case of GRB 071031, the C IV profile shows a smooth 
gradient with velocity once the —115 kms" 1 component has 
been subtracted off. 



ponents with W T > 200 mA occur at 1 000-5 000 kms -1 in 
20 per cent of GRB afterglow spectra. In our data, the in- 
cidence of HV C iv absorbers (6/7, or 86%) is much higher 
because (a) we have a higher sensitivity: the weakest HV 
Civ absorber in our sample shows W^ r =48 mA, and (b) we 
classify absorbers in the range 500-1000 kms" 1 as HV. 

The HV absorbers seen in GRB afterglow spectra do not 
represent a homogeneous population, as do the blueshiftcd 
absorption-line wings and (potentially) the strong, narrow 
absorbers at zgrb- On the contrary, we observe a broad 
range of properties among the HV absorbers. In three sight 
lines (toward GRBs 021004, 050730, and 060607), the HV 
components exhibit low-ionization absorption (clearest in 
C II A1334.532 and Sill A1260.422). However, the HV com- 
ponents toward GRBs 050820, 071031 and 080310 show no 
(or very weak) neutral-phase absorption, in either C n or 
Sill. So whereas the presence of neutral gas in the first 
three cases implies the absorbing material is not formed 
in the circumburst medium (because this medium is ex- 
pected to be fully ionize d by the burst; e.g.lVreeswiik et alj 
120071: IChen et all 120071 iProchaska et al.ll2008bft . we find 
that the properties of the HV components in the latter three 
cases are consistent with an origin in the close circumburst 
medium. This is particularly true for the HV gas toward 
GRB 071031, which is detected in the form of strong, O vi 
and C iv components that are not only absent in C II and 
Si II, but also show non-detections of Si iv. We have placed 
limits of iV(Civ)/V(Siiv) >110 at -560 kms" 1 , and >140 
at —510 kms -1 in the GRB 071031 spectrum. These ex- 
treme ratios are far higher than typical inter stellar values 
(Gala ctic ISM values for this ratio are «3; IZsargo et al.l 
120031 ) implying these HV components have a different, non- 
interstellar origin. 

Given the connection between GRBs and massive st ars 
(jMacFadven fc Wooslevl Il999l : IWoolsev fc Bloomll2006l) . a 
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leading explanation for the HV components in GRB after- 
glow spectra is that they trace outburst epi sodes from the 
mass ive-star progenitors of the GRBs (e.g. ISchaefer et al.l 
2003) , and thus arise i n the circumburst region. As pointed 
out by iBerger et alj (|2006h . a Wolf-Rayet outflow model 
would give rise to absorbers that are enriched in carbon and 
deficient in silicon, naturally explaining the high Civ/Siiv 
ratios measured in the HV absorbers. 

If the wind-blown bubbles around GRB progenitors are 
similar in size to those observed around Wolf-Rayet stars 
in the Milky Way, t heir typical radius will be «2-10 pc 
(|Gruendl et al.ll2000D . Several groups have recently mod- 
eled Wolf-Rayet winds and their rela tionship to the HV 
components i n GRB afterglow spect r a (iRamirez-Ruiz et alj 
120011 120051: IChevalier et all l2004t lEldridee et all 12006ft . 
and in particular, predictions for the high- ion signa- 
ture of Wol f -Raye t winds are given by Ivan Marie et al.l 
(120051 120071 120081 ). Depending on observing epoch and 
viewing angle, components at 150-700 kms -1 and 1000- 
2 000 kms " 1 are predicted at var ious stages of Wolf-Rayet 
evolution (jvan Marie et all 120051) : these velocities match 
the observed range of HV components in our sample. 
Recent models of the absorption-line signa tures of galac- 
tic wi nds (as opposed to sing le-star winds: Fangano et al.l 
120071 : iKawata fc RaucrJl2007t ISamui et al J 120081 1 generally 
predict high-ion components at velocities of only 100- 
200 kms -1 relative to the galaxy, though observations 
support the idea that galact ic winds can reach s peeds 
of »1000 kms" 1 J avejeLaLl I2002J; iTremonti et"aH 120071: 
Wci ner et al.ll2008l ). Nonetheless, the fairly narrow b- values 
measured in the HV components toward GRBs 071031 and 
080310 (not to mention the lack of Si iv) are at odds with 
the broader, c ollisionally ionized componen ts expected in 
galactic winds (|Qppenheimer fc Dav e 2006). Thus we con- 
clude that among the available models, Wolf-Rayet winds 
are the best explanation for the highly ionized HV compo- 
nents toward GRBs 071031 and 080310. 



5. Summary 

We have performed a systematic study of the high-ion ab- 
sorption near 2qrb in seven GRB afterglow spectra ob- 
served with UVES. Three of the seven were observed within 
15 minutes of the trigger by the Swift satellite. This search 
has provided high-resolution (6.0 kms -1 FWHM) profiles 
of absorption in the ions Ovi, Nv, Civ, Siiv, Siv, and 
S vi. The spectra contain a wealth of useful information 
on the properties of circumburst and interstellar gas in the 
GRB host galaxies. We have presented the results of Voigt 
profile fits to the high-ion absorption. Analysis of the line 
profiles shows that several types of high-ion absorption exist 
at zwzqrb, which we summarize in the following points. 

1. We detect strong high- ion absorption components ex- 
actly at .zgrb, always in O VI, C iv, and Si iv, usually in 
N v, and occasionally (where the data are unblended) in 
Siv and Svi. The high- ion column densities are signifi- 
cantly higher than those seen in the Milky Way ISM and 
in DLAs. We focus on N v since it is the least saturated 
and hence best-measured high ion. In three of seven 
afterglow spectra, the strong Nv absorption takes the 
form of a single component coincident in velocity with 
saturated absorption in the other hi gh ions, appear- 
ing to support the recent conclusion oflProchaska et all 



( 20083) that GRBs can photoionize N v in circumburst 
gas. However, in the remaining four cases, the Nv ab- 
sorption is either multi-component, offset from the other 
high-ions, or absent, and it is unclear whether the cir- 
cumburst hypothesis can explain these cases. In addi- 
tion, the Nv 6-values (median of 16 kms -1 ) are not 
narrow enough to rule out collisional ionization. Finally, 
for GRB 050730 we have placed a lower limit on the dis- 
tance to the S iv-absorbing gas of 400 pc, based on com- 
paring photo-excitation models to the observed column 
density of Siv and upper limit to Siv*. This supports 
an interstellar (rather than circumburst) origin for the 
strong high-ion absorption components at zgrb- 

2. In addition to the strong component at .zgrb, the Civ 
and Si iv absorption exhibits complex, multi-component 
profiles extending over several hundred kms -1 (exclud- 
ing the HV components, which we treat separately). 
The median total C iv line width in our seven spectra is 
280 kms -1 , similar to the value it takes in high-redshift 
DLA galaxies. This supports a galactic ISM origin for 
these multiple components. 

3. We detect asymmetric, blueshifted absorption-line 
wings in the Civ, Siiv, and (at lower S/N) Ovi pro- 
files in four of the seven GRB afterglow spectra. These 
wings are similar in shape and velocity extent (wlOO- 
150 kms -1 ) to positive-velocity wings seen in Ovi ab- 
sorption at zero redshift in the Milky Way halo, and two 
of the four wings show C iv/Si iv ratios of «3, equivalent 
to the value measured in the Milky Way ISM. We thus 
interpret the wings as tracing outflowing interstellar gas. 
The O vi column densities in the GRB wings are, in the 
mean, five times higher than the O vi column densities 
in the Milky Way wings (implying similarly higher mass 
flow rates), and the gradient dN a /dv is typically 5-10 
times higher. 

4. HV (500-5000 kms -1 ) components are seen in Civ in 
six of the seven spectra (all GRBs except 050922C). 
These components are not necessarily a single homo- 
geneous population, and many may arise in unrelated 
foreground galaxies. However, in the (previously unpub- 
lished) cases of GRBs 071031 and 080310, the ionization 
properties of the HV components (very high Civ/SilV 
ratios and an absence of neutral-phase absorption in Si 11 
or C 11) are consistent with a circumburst origin. Models 
of Wolf-Rayet winds from the GRB progenitor stars are 
capable of explaining the kinematics and ionization level 
of these HV components. 
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